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Remote  Observation  of  the  Spatial  Variability  of 
Surface  Waves  Interacting  With  an  Estuarine  Outflow 
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Abstract — This  paper  explores  the  application  of  phased -array 
high-frequency  (HF)  radars  to  identify  locations  of  enhanced  local 
waveheights.  Measurements  of  the  near-surface  current  velocities 
and  waveh eights  were  obtained  from  HF  radars  deployed  near  the 
mouth  of  the  Chesapeake  Bay  in  the  fall  ofl997.  The  radar-  derived 
near-surface  velocities  were  compared  with  the  upper  bin  (2-m 
depth)  of  four  upward -looking  acoustic  Doppler  current  profilers 
(ADCPs).  The  slopes  of  the  linear  correlations  w  ere  close  to  one  and 
the  rool-mean-square  (mis)  differences  were  similar  to  previous 
studies.  Significant  waveheight  {Hs)  estimates  from  both  radars 
were  compared  with  a  laser  height  gauge.  The  largest  differences 
were  observed  during  low  winds  due  to  overestimates  at  one  of 
the  radar  stations  and  during  storms  when  the  laser  measurement 
failed.  Further  analysis  focused  on  the  HF  radar  results  from  the 
more  reliable  of  the  two  sites.  The  mis  difference  between  this 
radar  and  the  in  situ  sensor  was  0,29  m.  Synoptic  observations  or 
II  #  over  the  Chesapeake  Bay  revealed  regions  of  current-induced 
wave  shoaling  and  refraction.  Ha  over  the  estuarine  outflow  in¬ 
creased  between  19-50%  relative  to  the  incident  Hh  in  light  on¬ 
shore  winds  (~5  m/s).  In  stronger  winds  (>I0  m/s)*  H #  also  in¬ 
creased  by  up  to  25%  when  there  was  a  tidal  outflow  in  the  surface 
layer*  although  the  near-surface  currents  were  responding  to  both 
the  wind  and  the  ebbing  tide.  H *  was  not  enhanced  when  the  out¬ 
flow  was  below  a  thicker  layer  (>  5  ml  of  wind -forced  onshore  flow  . 

Index  Terms — Coastal  currents*  high-frequency  (HF)  radar* 
shoaling*  wave  refraction. 


L  Introduction 

AS  MANY  sailors  leaving  a  coastal  inlet  during  an  ebb  tide 
have  observed,  the  local  waveheight  often  increases  dra¬ 
matically  where  the  outgoing  flow  meets  incoming  waves.  Estu¬ 
arine  outflows  often  exhibit  regions  of  significant  current  shear* 
the  location  of  which  depends  on  discharge,  tidal  stage*  and 
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wind  forcing.  Surface  waves  propagating  over  these  variable 
currents  and  topography  experience  both  refraction  and  shoaling 
which  can  lead  to  spatially  inhomogeneous  waveheights. 

Phased-array  high-frequency  (HF)  radar  systems  offer  the 
promise  of  making  synoptic  measurements  of  both  currents  and 
waveheights  simultaneously  over  a  coastal  region  of  interest, 
therefore  allowing  the  near-real-time  identification  of  regions 
where  enhanced  waves  pose  a  hazard  to  maritime  activities. 
Furthermore,  the  high  resolution  1  km)  surface  current  and 
wave  information  can  provide  the  necessary  information  to  in¬ 
corporate  into  numerical  models  of  these  processes.  The  focus 
of  this  paper  is  on  the  first  of  these  capabilities*  the  observation 
of  regions  of  locally  increased  waveheights  due  to  shoaling  and 
refraction  over  topography  and  highly  sheared  currents. 

Section  II  of  this  paper  provides  background  information 
on  wave  refraction  and  shoaling  over  variable  currents.  HF 
radar  observations  of  waves  and  currents  are  then  discussed. 
The  radar  measurements  of  both  currents  and  waves  during  the 
third  Chesapeake  Outflow  Plume  Experiment  (COPE-3)  are 
compared  with  available  in  situ  observations  to  establish  their 
accuracy.  Two  cases  when  the  wind  direction  was  opposing 
the  ebb-tidal  outflow  from  the  Chesapeake  Bay  are  examined 
in  Sections  III  and  IV.  The  first  of  these  was  a  case  with  light 
winds  (^6  m/s)*  while  the  second  had  winds  >12  m/s.  The 
source  of  the  observed  spatial  variability  of  the  wave  field  in 
these  two  situations  is  explored  in  Section  IV. 

A,  Wave  Refraction  and  Shoaling 

Surface  waves  are  affected  by  many  processes  as  they  prop¬ 
agate  across  a  shallow  shelf.  In  particular*  wave  energy  is 
dissipated  by  wave  breaking  and  bottom  friction  [1].  The  local 
waveheight  is  increased  by  shoaling  caused  by  decreasing  water 
depths.  Waveheights  also  increase  when  waves  are  moving 
against  an  opposing  current  and  decrease  in  a  following  current. 
Wave  shoaling  in  opposing  currents  depends  on  their  wave  length 
and  direction  relative  to  the  current  direction  and  magnitude. 
Focusing  or  divergence  of  wave  rays  by  refraction  can  also  be 
caused  by  either  depth  or  current  gradients.  Because  of  the  many 
processes  involved*  there  can  be  significant  spatial  and  temporal 
variability  of  waveheights  in  tidally  dominated  coastal  regions. 

If  dissipative  processes  are  neglected*  the  shoaling  of  nor¬ 
mally  incident  surface  gravity  waves  propagating  over  an  op¬ 
posing  current  and  a  sloping  beach  can  be  estimated  using  linear 
theory  by  H2/Hx  =  yf{C9 1  -  V\)/{Cg2  -  U2).  where  H\  and 
H2  are  the  waveheights*  C9\  and  Cy2  are  the  wave  group  veloc¬ 
ities*  and  U\  and  U2  are  the  currents  in  the  wave  direction  at 
two  points  along  the  wave  propagation  path  [2J.  From  this  ex¬ 
pression*  it  can  be  seen  that  topography  and  currents  influence 
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Fig,  1 .  (■)  Critical  velocity  at  which  surface  waves  grow  without  bound  related  to  the  deep-w  ater  wave  frequency.  (* )  Fraction  of  energy  of  fetch  limited  waves 
with  wind  of  10  m/s  and  inverse  wave  age  of  1 ,673  that  occur  at  higher  frequencies  than  shown  [36J.  (-)  Energy  fraction  of  fully  developed  sea  ai  10  m/s  that  Occurs 
at  higher  frequencies. 


the  shoaling  of  surface  waves  in  distinctly  different  ways.  Small 
wave  number  (fc)  waves  are  most  affected  by  topography  be¬ 
cause  they  “feel”  the  bottom  before  higher  wave  number  waves. 
Conversely,  the  high  wave  number  portion  of  the  gravity  wave 
spectrum  is  most  strongly  impacted  by  currents  because  the 
waves  have  lower  group  velocities. 

The  forward  propagation  of  short  gravity  waves  (k  >6)  will 
he  completely  halted  by  currents  greater  than  I  m/s.  This  critical 
opposing  current  velocity  { Ur )  for  a  particular  k  is  given  by 
Uc  —  Cg a/2,  where  C3 o  is  the  deep-water  incident  wave  group 
velocity  [3].  For  fully  developed  seas  under  moderate  winds, 
much  of  the  wave  energy  will  be  at  wave  numbers  such  that  the 
wave  energy  will  not  approach  Uc  for  currents  less  than  1,5  m/s 
(Fig.  1).  However,  for  light  winds,  young  seas,  or  fetch  limited 
conditions,  a  significant  proportion  of  the  wave  energy  occurs  at 
frequencies  where  unbounded  growth  will  occur  in  an  opposing 
current  less  than  1.5  m/s. 

Because  the  speed  of  wave  propagation  changes  as  waves 
move  across  variable  currents  or  topography,  the  direction 
of  propagation  of  the  waves  will  also  change  by  Snell's  Law 
Sin  (#i)/{Ci  —  t/i)  -  Sin  (^)/(Ca  -  TV)*  For  a  uniformly 
sloping  beach,  this  relationship  causes  waves  incident  at  an 
angle  to  the  beach  to  turn  toward  the  shore- normal  as  the  depth 
decreases  and  the  waves  slow  down.  For  waves  propagating 
obliquely  into  a  sheared  opposing  current  the  same  effect 
will  occur.  The  wave  rays  will  bend  toward  the  normal  to 
the  streamlines  in  the  direction  of  increasing  current  velocity. 
Kenyon  [4]  demonstrated  that  in  cases  of  weakly  sheared  flow 
the  ray  curvature  is  equal  to  the  ratio  between  the  flow  vor- 
ticity  and  the  wave  group  velocity.  For  real  coastal  conditions 
with  complicated  topography  and  current  patterns,  the  local 
waveheights  can  be  significantly  affected  by  the  focusing  and 
defocusing  of  the  wave  energy  through  refraction. 


The  equations  governing  the  refraction,  diffraction,  shoaling, 
and  dissipation  of  wave  energy  through  bottom  friction  and 
breaking  are  highly  nonlinear.  In  general,  approximations  to  the 
full  equations  are  used  to  estimate  local  wave  conditions  given 
known  incident  waves.  Kirby  and  Dairy mple  [5]  developed  a 
numerical  solution  to  a  parabolic  approximation  of  the  equa¬ 
tions,  which  has  been  widely  employed  to  model  coastal  wave 
conditions.  Numerical  refraeti on/d iffracl ion  models  have  heen 
well  tested  in  the  field  for  relatively  simple  topographies  and 
over  shoals  and  barred  beaches  in  laboratories  [6],  However, 
providing  sufficient  spatial  distribution  of  measurements  for 
quantitative  calibration  and  validation  of  numerical  models  is 
difficult  for  many  coastal  applications. 

B ,  HF  Radar  Measure  me  n  rs 

HF  radars  can  sample  surface  currents  over  large  enough 
areas  with  sufficient  spatial  resolution  to  resolve  tidal  circula¬ 
tion  patterns  and  estuarine  outflow  plumes.  Phased-array  HF 
radars  also  can  observe  the  shoaling  and  refraction  of  surface 
waves  over  the  same  areas.  The  combination  of  these  wave  and 
current  observations  offers  opportunities  for  real-time  mapping 
for  maritime  uses  as  well  as  assimilation  into  or  validation  of 
numerical  circulation  and  wave  models. 

Most  HF  radars  operate  in  a  monostatic  mode  of  operation,  in 
that  they  transmit  a  pulsed  or  continuous  wave  signal  and  then 
receive  the  back  scattered  signal  at  the  same  location.  The  time 
required  for  the  signal  to  travel  a  particular  distance  over  the 
ocean  and  return  to  the  receiver  can  be  precisely  determined  and, 
therefore,  localizing  the  observations  in  range  is  easily  done. 

To  localize  the  received  signal  in  azimuth  is  a  more  com¬ 
plicated  matter.  Phased-array  HF  radar  systems  observe  the 
in-phase  (7)  and  quadrature  (Q)  signals  at  multiple  antennas 
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Doppler  Frequency  [Hz] 

Fig.  2.  Sample  echo  Doppler  spectrum  derived  from  the  ocean  surface  current  radar  (OSCR)  HF  radar  operating  ai  25.4  MHz.  First-order  Bragg  peaks  are  shown. 
The  range  of  frequencies  over  which  the  second-order  energy  is  integrated  to  derive  waveheight  is  denoted  by  low  and  high. 


linearly  spaced  at  half  the  radar  wavelength  and  apply  prede¬ 
termined  phase  relationships  to  obtain  azimuthal  localization. 
The  theoretical  resolution  of  the  azimuth  angle  (0)  is  equal  to 
the  ratio  between  the  radar  wavelength  (A)  and  the  length  of 
the  receiver  array  (D).  For  the  phased-anray  system  used  in  this 
paper 

(5) = (5$  -  *  - r-6’' 

The  back  scattered  signal  obtained  from  this  beam  forming  tech¬ 
nique  is  recorded  for  the  duration  of  the  transmission  at  each 
range  bin.  Enough  observations  of  the  backscattered  signal  from 
the  ocean  surface  are  obtained  during  each  sample  period  to 
calculate  the  echo  Doppler  spectrum  (Fig,  2)  from  the  Fourier 
transform  of  the  combined  /  and  Q  signals. 

As  the  number  of  observing  antennas  is  increased,  the 
signal-to-noise  ratio  (SNR)  of  the  Bragg  scattering  signal  is 
also  increased.  This  renders  the  positions  of  the  first-order 
Bragg  peaks  relatively  insensitive  to  propagation  path  distor¬ 
tions  of  the  signal  for  a  16-element  receive  array  as  used  here. 
However,  under  some  circumstances  such  as  when  the  receive 
array  is  placed  at  a  large  angle  to  the  shoreline  the  antenna 
beam  pattern  can  become  asymmetrical  [7],  This  has  been  ob¬ 
served  to  corrupt  the  velocity  signal  when  the  beam  distortion 
was  combined  with  a  significant  spatial  in-homogeneity  of  the 
radar  cross  section  of  the  ocean  surface  [8),  Most  deployments 
of  phased-array  radars  have  not  encountered  such  conditions 
and  have  not  been  significantly  impacted  by  beam  pattern 
distortions. 

Direction  finding  techniques  use  compact  antennas  and  do  not 
construct  the  echo- Doppler  spectrum  at  each  range  and  azimuth 
bin.  Rather  they  identify  the  direction  that  signals  at  specified 


frequencies  are  emanating  from.  The  azimuthal  resolution  of 
HF  radars  using  direction  finding  is  limited  by  the  short  aper¬ 
ture  length  and  they  are  more  sensitive  to  antenna  beam  pattern 
distortions.  In  many  cases,  the  actual  beam  pattern  must  be  mea¬ 
sured  using  boat  mounted  transponders  to  obtain  optimal  results 
with  direction  finding  systems.  The  echo-Doppler  spectrum  can 
be  constructed  to  make  wave  observations  |9],  but  this  typically 
requires  averaging  over  large  azimuth  angles, 

i)  Near-Surface  Currents:  HF  radar  signals  are  preferen¬ 
tially  scattered  off  waves  with  lengths  of  half  the  transmitted 
wavelength  by  the  Bragg  scattering  mechanism  [lOj.  This  pro¬ 
duces  a  dominant  first-order  peak  in  the  echo-Doppler  spectrum 
(Fig,  2).  As  the  surface  waves  that  are  responsible  for  the  Bragg 
scattering  peaks  are  advected  either  toward  or  away  from  the 
radar  by  near-surface  currents  their  relative  speed  is  shifted. 
The  surface  current  velocity  is  then  determined  by  the  offset 
(Doppler  shift)  of  the  first-order  peaks  from  the  zero  current 
condition.  The  resolution  of  the  radial  velocity  measurements  is 
determined  by  the  precision  with  which  the  frequency  shift  can 
be  determined  from  the  Doppler  spectrum.  For  a  pulsed  system, 
the  frequency  resolution  is  dependent  upon  the  pulse  length  and 
the  duration  of  sampling  which  determine  the  total  number  of 
pulses  used  to  derive  the  Doppler  spectrum. 

When  the  radial  measurements  from  two  independent  stations 
are  combined  to  estimate  the  vector  current,  the  measurement 
accuracy  is  reduced  when  the  radialsare  not  orthogonal.  The  op¬ 
timal  vector  resolution  is,  therefore,  dependent  upon  the  angle 
between  the  observations  and  the  true  current  direction  [II], 
[12].  The  realized  accuracy  is  dependent  upon  atmospheric  con¬ 
ditions,  sea  state,  and  the  signal  properties  observed  at  each  of 
the  receiver  antennas. 

HF  radars  have  resolved  the  highly  sheared  current  fields 
associated  with  the  inshore  edge  of  the  Florida  current  [13], 


m 

small-scale  instability  vortices  that  occur  in  the  coastal  buoy- 
ancy  current  along  the  Outer  Banks  of  North  Carolina  f  14],  and 
the  tidal  circulation  patterns  associated  with  the  Chesapeake 
Bay  outflow  plume  [15],  [16]. 

2)  Surface  Wave  Parameters:  The  echo-Doppler  spectrum 
contains  information  on  both  the  surface  currents  and  the  sur¬ 
face  waves.  In  addition  to  the  first-order  Bragg  resonances,  a 
typical  spectra  contains  significant  energy  at  frequencies  around 
the  first -order  peaks  (Fig.  2).  These  second-order  returns  contain 
backscattered  energy  resulting  from  multiple  reflections  of  the 
radar  signal  as  well  as  the  hydrodynamic  combination  of  surface 
waves  to  produce  Bragg  scattering  [17].  This  second-order  en¬ 
ergy,  therefore,  contains  information  on  the  surface  wave  spec¬ 
trum,  This  information  is  ideally  available  at  every  range  and 
azimuth  bin  over  the  sampling  region.  However,  the  range  of 
useful  wave  information  is  more  limited  than  for  current  mea¬ 
surement  because  of  the  lower  SNR  associated  with  the  second- 
order  peaks.  Phased-array  HF  radars  can  provide  sufficient  SNR 
and  azimuthal  resolution  to  allow  spatial  localization  of  wave 
observations. 

Methods  to  invert  the  second-order  contribution  to  the 
Doppler  spectrum  to  estimate  the  surface  wave  directional 
spectrum  have  been  developed  by  [18]  and  [19],  In  [20],  it  was 
demonstrated  that  inversion  methods  can  obtain  reliable  wave 
estimates  when  compared  with  moored  buoys.  These  methods 
require  that  either  strict  bounds  are  placed  on  the  frequency 
range  of  the  inversion  [18]  or  that  the  raw  quadrature  signal  are 
collected  and  stored  [19]. 

Historically,  for  field  operations  the  large  amount  of  disk 
storage  required  to  archive  the  raw  data  was  often  difficult  to 
manage.  Archiving  the  Doppler  spectra  derived  from  the  1  and 
Q  signals  using  fast  Fourier  transforms  required  a  factor  of 
three  less  storage  and  was  often  the  only  reasonable  alternative 
for  waveheight  observation.  Modem  computational  resources 
have  effectively  eliminated  these  restrictions,  however,  statis¬ 
tical  parameters  of  the  wave  field  such  as  Hs  are  still  widely 
used  to  characterize  wave  conditions. 

Empirically  based  methods  to  derive  omnidirectional  wave- 
height  estimates  from  the  second-order  energy  were  originally 
developed  by  [17],  This  type  of  approach  has  been  tested  and 
validated  for  the  phased-array  OSCR  HF  radars  by  [21  ]  and  has 
been  proven  to  provide  reliable  estimates  of  Hs*  The  use  of  in¬ 
formation  from  overlapping  stations  provides  two  estimates  of 
Hs  at  each  measurement  location.  The  independent  measure¬ 
ments  can  then  be  compared  to  estimate  the  reliability  of  the 
observations.  No  directional  information  can  be  obtained  from 
a  single  radar  alone  without  averaging  over  large  azimuth  an¬ 
gles,  however  by  blending  radar  observations  and  models  [22] 
estimated  directional  spectra  from  a  single  site. 

The  simple  empirical  approach  was  employed  here  to  provide 
synoptic  mapping  of  Hs.  The  waves  and  currents  were  simul¬ 
taneously  observed  with  the  OSCR  phased-array  radar  and  in 
situ  sensors  during  COPE-3.  The  accuracy  of  the  remote  obser¬ 
vations  as  established  by  the  in  situ  sensors  is  discussed  in  this 
paper  in  detail.  The  spatial  variability  of  the  near-surface  cur¬ 
rents  and  waveheights  and  the  interaction  between  them  during 
conditions  when  there  was  significant  horizontal  and  vertical 
shear  is  then  explored. 
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II.  COPE-3  Experiment 

The  OSCR  system  mapped  surface  currents  and  waveheights 
over  the  inner-shelf  offshore  the  Chesapeake  Bay  mouth  during 
the  COPE-3  experiment  in  fall  1997  (Fig,  3).  Vector  current 
fields  were  obtained  for  the  period  from  yearday  (YD)  287  to 
YD  334.  One  station  (HF-N)  was  located  at  the  Fort  Story,  VA, 
U.S.  Army  base  just  to  the  south  of  the  estuary  mouth.  The 
second  radar  station  (HF-S)  was  located  18  km  to  the  south¬ 
east  at  the  U.S.  Navy  Fleet  Combat  Training  Center  Atlantic  in 
Dam  Neck,  VA, 

Both  stations  were  operated  at  25.4  MHz  and  transmitted  a 
pulsed  signal  with  a  wavelength  of  1 1.6  m.  The  corresponding 
ocean  wave  Bragg  scattering  wavelength  was  5,8  m.  This  pro¬ 
vided  an  estimate  of  the  bulk-average  of  the  currents  over  the 
top  0.4-1 .5  m  of  the  water  column.  The  horizontal  bin  size  was 

1  km2  over  an  area  of  approximately  30  km  x  44  km  (Fig,  3), 
There  was  a  5-min  period  of  transmission  at  HF-S  followed  by 
5  min  at  HF-N  and  then  10  min  for  processing  the  received  sig¬ 
nals.  This  sampling  provided  468  pulses  to  produce  the  echo- 
Doppler  spectrum  from  which  the  radial  velocities  were  deter¬ 
mined  with  an  optimal  accuracy  of  0.02  m/s. 

Estimates  of  Hs  were  obtained  at  each  measurement  location 
through  the  ratio  of  the  second-order  backscattered  energy  to  the 
first-order  Bragg  scattering  energy  following  [17].  This  empir¬ 
ical  approach  used  linear  scaling  coefficients  to  minimize  the 
mean  square  difference  between  HF  radar  observed  and  in  situ 
Hs  measurements  for  Hs> 0,53  m.  The  coefficients  used  to 
scale  the  second-order  spectrum  for  this  paper  were  calibrated 
from  earlier  OSCR  measurements  at  Duck,  NC  [21  ]. 

Data  from  five  acoustic  Doppler  current  profilers  (ADCP) 
were  used  in  the  validation  studies  (Fig.  3).  The  ADCPs  were 
307.2-kHz  broadband  * ‘workhorse  sentinels"  manufactured  by 
RD  Instruments,  San  Diego,  CA.  They  were  deployed  in  bottom 
mounts  in  upward-looking  configurations  on  October  6,  1997 
(YD  279)  and  recovered  on  November  10,  1997  (YD  314).  The 
vertical  bin  size  was  I  ni  and  valid  currents  were  acquired  from 

2  m  beneath  the  surface  to  about  4  m  above  the  bottom.  Water 
depths  at  the  moorings  ranged  from  10  to  18  m.  The  sampling 
scheme  used  was  a  1-min  burst  of  120  samples,  which  was  re¬ 
peated  at  5-min  intervals.  Burst  sampling  was  used  to  mini¬ 
mize  aliasing  of  the  measurement  by  surface  gravity  waves.  Ve¬ 
locity  data  from  each  sample  were  resolved  into  northward  and 
eastward  components,  then  averaged  over  the  burst.  Raw  data 
were  further  smoothed  with  a  three- point  ru nning  mean  before 
analysis. 

Wind  and  in  situ  wave  measurements  were  obtained  from 
the  National  Oceanic  and  Atmospheric  Administration  (NOA A) 
Coastal  Marine  Automated  Network  (CMAN)  Chesapeake  Bay 
Lighthouse  Tower  (CHLV2)  station.  The  station  was  located 
within  the  HF  radar  measurement  domain  at  a  distance  of  25  km 
from  HF-N  and  26  km  from  HF-S  (Fig,  3),  The  anemometer 
was  located  at  the  top  of  the  tower  at  a  height  of  43.3  m  above 
the  mean  water  level.  Wind  measurements  were  converted  to 
standard  10- m  height  by  assuming  a  log  profile.  The  wind  ve¬ 
locity  and  waveheight  were  recorded  every  10  min  throughout 
the  experiment.  The  waveheight  was  derived  from  a  laser  height 
gauge  mounted  on  the  tower  [23]  and  was  nondirectional 
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Fig.  3.  Overview  of  Mid-Atlantic  region  (right).  Topography  and  measurement  locations  for  the  COPE-3  experiment  shown  as  inset.  Depth  contours  are  shown 
every  5  m  The  center  of  each  HF  radar  sample  region  is  shown  as  (.),  The  shore  stations  are  denoted  by  large  solid  circles.  The  ADCP  moorings  AL  AZ  A3.  A4. 
and  A5  are  shown  with  triangles.  The  CHLV2  CM  AN  station  is  shown  as  a  square. 


TABLE  1 

Comparisons  of  OSCR  HF  Radar  Current  Velocities  With  2-m  Bin  of  ADCP,  a  Is  the  Standard  Deviation  of  the  HF  Radar  Component.  AU  and 
AV  ARE  THE  DIFFERENCES  BETWEEN  EAST-WEST  AND  NORTH-SOUTH  COMPONENTS,  RESPECTIVELY.  CORRELATION  SLOPE (m)  AND  BIAS  (b)  FOR  EACH 
Component  as  Well  as  the  Complex  Correlation  Coefficient  (p)  and  Phase  <^)  Are  Shown 


U 

V 

m/s 

rms 

U 

m/s 

m 

b 

m/s 

m/s 

rms  V 

m/s 

m 

b 

m/s 

A2 

0.25 

0.11 

LI 

0.05 

0.31 

0.10 

0.9 

*0.02 

0,94 

-9 

A3 

0.22 

0.09 

1.0 

0.02 

0.21 

0,09 

LI 

-0.02 

0.91 

-4 

A4 

0,18 

0,07 

1,0 

0.02 

0.19 

0,08 

L2 

-0.01 

0.90 

2 

A5 

0.16 

0,07 

LI 

-0.0 1 

0.21 

0.10 

L2 

-0.02 

0.94 

12 

in.  Results 

A ,  Surface  Current  Comparisons  With  ADCPs 

The  East- West  (U)  and  North-South  (V)  velocity  compo¬ 
nents  of  each  ADCP  were  compared  with  the  nearest  HF  radar 
measurement.  The  HF  radar  U  velocities  at  A 1  had  unaccept¬ 
ably  large  errors  because  of  the  targe  angle  between  the  two  ra- 
dials  [11].  The  results  from  this  location  were,  therefore,  omitted 
from  the  analysis.  The  second  bin  from  the  surface  of  the  ADCP 
was  used  in  each  case  because  of  the  contamination  of  the  first 
bin  due  to  sidetobe  reflections. 


The  statistics  of  the  comparisons  between  the  unfiltered 
OSCR  observations  in  the  near-surface  layer  and  the  ADCP 
revealed  good  agreement  between  the  sensors  (Table  I),  The 
slopes  of  the  best  fit  regression  curve  were  either  1,1  or  1  in 
each  case  for  U  with  biases  below  0.05  m/s  [24].  The  slopes 
of  the  V"  comparison  ranged  from  0.9  to  1 .2  with  biases  be¬ 
tween  -0.009  m/s  and  -0,018  m/s.  The  complex  correlation 
coefficient  was  0.9  or  above  at  each  of  the  four  moorings. 

The  rms  difference  between  the  HF  radar  and  the  ADCP 
U  component  at  A2  was  larger  than  at  the  other  moorings 
(0.1 1  m/s).  This  was  in  a  region  with  large  tidal  excursions  and 
this  difference  as  a  fraction  of  the  velocity  range  was  about 
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Fi  g .  4  .  Ti  me  series  o  f  u  n  fi  ltered  hou  rly  su  rface  ( •  ,  OSCR )  and  A  DC  P  obse  rved 
subsurface  (solid)  velocities  ai  the  A2  mooring,  (a)  T  (Easi-West)  component 
positive  toward  cast,  (b)  V  (North-South)  component  positive  toward  north. 
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Fig.  5.  Time  series  of  the  significant  waveheighi  observed  by  HK  radars  HF-S 
(X).  HF^N  (+).  and  CHLV2  (*  J  tower  mounted  laser  during  the  COPE  3  ex  per 
iment  in  1997. 


1 1  %.  The  time-series  comparison  of  these  velocities  reveals  the 
dose  relationship  between  the  two  measurements  (Fig.  4).  The 
only  observable  difference  is  that  the  OSCR  U  observations 
were  consistently  larger  than  the  ADCP  values  at  the  peak  of 
the  ebb  tide  and  were  intermittently  larger  at  the  peak  of  the 
flood  tide.  The  amount  of  difference  varied  over  each  tidal 
cycle  but  was  typically  between  0.05  and  0. 1 5  m/s. 

B.  Comparisons  of  HF  Radar  Waveheighi s  With  Laser  Height 
Gat4ge 

The  H s  observed  by  HF-S  and  HF-N  at  the  closest  measure¬ 
ment  cell  were  compared  with  those  obtained  at  CHLV2.  No 
adjustments  to  the  scaling  coefficients  obtained  from  [21]  were 
made  to  the  values  used  here.  The  radar  measurements  track 
the  passage  of  mid  latitude  fronts  in  a  similar  manner  to  CHLV2 
(Fig.  5k  although  CHLV2  experienced  data  loss  at  the  peak  of 
the  most  energetic  storms. 

The  largest  waveheights  during  the  experiment  were  recorded 
over  the  period  from  YD  291  to  YD  295.  The  only  measurement 
collected  throughout  this  period  was  from  HF-S,  HF-S  closely 
matched  HF-N  during  the  period  before  the  peak  and  HF-S 
matched  CHLV2  as  the  waveheighi  decreased. 

On  YD  312,  both  of  the  radars  and  CHLV2  recorded  Ha 
values  between  2.5-2, 8  m  with  the  HF-S  reading  slightly  higher 
than  the  other  two.  Later  on  YD  317,  all  three  measurements 
again  tracked  the  passage  of  a  storm  system  with  Hs  exceeding 
3  m.  In  this  case,  both  radars  recorded  peak  values  larger  than 
CHLV2,  hut  again  there  were  dropouts  in  the  laser  measure¬ 
ments  at  the  peak  of  the  storm. 

One-sided  regressions  revealed  an  rms  difference  of  H  s  be¬ 
tween  HF-S  and  CHLV2  of  0,29  m  (Table  II),  The  rms  dif¬ 
ference  between  HF-N  and  CHLV2  was  considerably  larger 
(0.38  m)  and  it  was  0.41  m  between  the  two  radar  sites.  The 
two-sided  regressions  produced  smaller  rms  differences  from 
0.09  to  0. 14  m  in  each  case  (Table  II).  Inspection  of  the  scatter 
diagrams  (Fig,  6)  revealed  that  the  main  discrepancies  were 
large  values  of  Hs  recorded  at  HF-N  when  there  were  only 


small  waves  observed  in  the  other  two  measurements  [Fig.  6(h) 
and  (c)j.  The  HF-N  intermittently  recorded  Hs  values  0.5-1  m 
larger  than  HF-S  and  CHLV2  during  periods  with  otherwise 
low  waveheights.  The  large  error  at  small  waveheights  in  HF-N 
is  reflected  in  percent  relative  error  values  exceeding  10%  for 
waveheights  less  than  t  m  [Fig.  6(d)],  The  HF-S  percent  rela¬ 
tive  error  did  not  exceed  10%  except  for  waveheights  less  than 
0.3  m  [Fig.  6(d)], 

The  intermittent  large  HF-N  observed  Hs  values  occurred 
most  often  between  YD  303— YD  310  and  from  YD  320  to  the 
end  of  the  experiment.  During  the  last  period,  the  peaks  occurred 
almost  daily,  but  there  was  no  consistent  periodicity  to  their  oc¬ 
currence.  The  events  of  large  Hs  at  HF-N  typically  persisted  for 
2—4  h.  The  echo- Doppler  spectrum  recorded  during  these  events 
contained  large  amounts  of  energy  far  removed  from  the  Bragg 
peaks  (Fig.  7).  This  indicated  that  outside  signals  contaminated 
the  observations  at  these  times.  The  HF-S  observations  did  not 
have  these  short  duration  peaks  except  for  a  single  event  on  YD 
304. 

The  hourly  radar  wave  observations  from  HF-S  (767  9(H) 
points)  were  compared  with  the  simultaneous  current  measure¬ 
ments  at  the  same  location.  There  was  no  correlation  between 
the  velocity  magnitude  (Vinag)  or  either  of  the  components  and 
the  observed  Hs  (r  —  -0.03*  -0. 10,  and  -0,20  for  Vinag,  IF 
and  V).  Similar  correlation  values  were  obtained  with  smaller 
subsets  located  near  the  CHLV2  tower  and  with  the  laser  gauge. 
The  slightly  higher  correlation  with  the  V  component  (-0,2) 
was  related  to  the  physical  link  between  wind/wave  forcing  and 
this  component  of  the  currents. 

The  lack  of  correlation  between  the  radar  observations  of  Hs 
and  cunem  velocity  for  the  full  data  set  demonstrated  that  the 
observed  wave -current  interaction  during  particular  events  was 
not  an  artifact  of  the  radar  measurement.  This  was  expected 
because  the  two  measurements  were  derived  from  the  echo- 
Doppler  spectrum  in  fundamentally  distinct  ways.  The  surface 
current  was  proportional  to  the  displacement  of  the  first-order 
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TABLE  II 

One-Sided  Regressions  of  H *  Measurements  From  OSCR  HF  Radar  With  CHLV2  CM  AN  Observations 


rms  diff.  Hs  (m) 

Slope  of  best  fit 

Bias  of  best  fit  (m) 

HF-N  versus  CHLV2 

0.38 

0.77 

0.31 
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Fig.  6.  fa)  Significant  waveheight  as  observed  from  the  CHLV2  CM  AN  station  versus  the  closest  OSCR  HF  radar  cell  (# )  as  observed  from  HF-S  for  ihc  duration 
of  COPE-3,  Solid  line  is  the  unit  slope  line.  Dashed  line  is  the  best  fit  two-sided  regression,  (b)  Significant  wavcheight  as  observed  from  CHLV2  versus  the  closest 
OSCR  HF  radar  cell  (#)  as  observed  from  HF-N  for  the  duration  of  COPE-3.  Solid  line  is  the  unit  slope  line.  Dashed  line  is  the  slope  of  the  best  fit  two-sided 
regression,  te)  Ms  as  observed  from  the  HF-S  versus  the  HF-N  (f).  Solid  line  is  the  unit  slope  line,  dashed  line  is  the  two-sided  regression  best  fit,  (d)  Relative 
error  percentage  as  a  function  of  Hs,  Solid  line  is  CHLV2  versus  HF-S  Dashed  line  isCHLV2  versus  HF-N.  Dotted  line  is  HF-S  versus  HF-N. 
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Fig.  7,  Examples  of  large  spurious  signal  in  OSCR  echo- Doppler  spectra  observed  at  HF-N  during  COPE-3  (a)  YD  292.87  and  (b)  YD  290.96 
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Fig,  8.  Lunar  semidiurnal  (M2)  tidal  constituent  ellipses  derived  from  har¬ 
monic  analysis  of  45-d  records  of  OSCR  vectors  during  COPE-3.  Horizontal 
axis— longitude  East  (decimal  degrees).  Vertical  axis — latitude  North  (decimal 
degrees). 

peaks  in  frequency  space,  while  the  Hs  was  proportional  to  the 
ratio  of  the  backscattered  power  in  the  second -order  region  to 
the  first -order  peak  (Fig.  2). 

C  Tula!  Currents 

Analysis  of  the  longest  continuous  segment  of  the  CGPE-3 
surface  current  observations  (35  d)  confirmed  that  the  lunar 
semidiurnal  (M2)  was  the  most  energetic  tidal  component.  There 
was  little  vertical  structure  of  the  M2  constituent  (not  shown) 
confirming  that  it  was  primarily  barotropic  1 16],  The  major  axes 
of  the  M2  ellipses  were  generally  aligned  with  the  Chesapeake 
Bay  mouth  (Fig.  8),  The  strength  of  the  M2  tide  decreased  with 
distance  from  the  mouth.  This  decreasing  velocity  of  the  tidal 
flow  with  distance  from  the  estuary  mouth  caused  surface  current 
convergences  and  divergences  on  particular  phases  of  the  tide. 

During  the  ebb  tide,  the  outflow  of  estuarine  water  from  the 
Chesapeake  Bay  typically  turned  to  the  southeast  and  moved 
through  the  measurement  domain.  The  highest  surface  velocities 
were  usually  associated  with  the  ^20-m  deep  shipping  channel 
in  the  northwest  comer  of  the  domain  (Fig,  3).  The  buoyant 
plume  associated  with  the  estuarine  outflow  w  often  formed  a 
surface  convergence  front  where  it  overrode  denser  shelf  water. 

D.  Waves  Propagating  Into  Opposing  Currents 

As  the  ebb  tide  exited  the  bay  mouth  on  YD  287.9,  the  wind 
was  blowing  from  the  southeast  at  5,8  m/s.  Under  the  assump¬ 
tion  that  the  waves  were  propagating  in  the  wind  direction,  this 
produced  a  situation  in  which  the  locally  generated  wind  waves 
encountered  a  distinct  surface  current  gradient  (Fig.  9).  Offshore 
the  currents  were  weak  (^0.15  m/s)  and  rotated  to  the  right  of 
the  wind  direction.  In  the  region  of  the  estuary  mouth,  the  sur¬ 
face  currents  reached  magnitudes  over  1  m/s  in  opposition  to  the 
wind  and  wave  direction. 
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Fig,  0,  Surface  currem  (armws.  meter  per  second  scale  shown  in  lower  right 
comer)  vectors  from  combination  of  HF-S  and  HF-N  radial*  and  Ha  (contours, 
meter)  as  observed  by  HF-S  on  YD  287.78.  Wind  velocity  at  CHLV2  is  shown 
in  the  upper  left.  Only  those  values  obtained  ill  ranges  <=  2b  kin  from  HF-S 
are  shown. 


Fi g,  1 0.  Freq ue nc y  s pectra  observed  at  C H LV 2 .  (a )  Y I>  287 .  7  8 .  ( b }  Y  D  305 . 33. 


Over  the  southern  half  of  the  measurement  domain,  the 
HF-S-observed  Hs  was  generally  <0,4  m  (Fig,  9).  Hs  increased 
to  0.6-0. 8  m  in  the  region  of  opposing  currents  near  the  A 2 
mooring.  The  largest  Hs  values  were  associated  with  the  strong 
outflow  in  the  shipping  channel  that  was  directly  in  opposition  to 
the  wind/wave  direction.  In  the  center  of  the  domain  where  the 
surface  currents  were  directed  more  offshore,  the  Hs  values  were 
less  than  0.4  m.  The  frequency  spectrum  recorded  at  CHLV2  at 
this  time  revealed  that  in  addition  to  the  wind  wave  contribution 
around  0.17  Hz,  there  was  a  strong  swell  peak  in  the  spectrum 
at  0.07  Hz  Fig.  10(a).  Inspection  of  the  directional  spectrum 
recorded  11 1  km  east  of  the  measurement  region  revealed  that 
the  swell  was  incident  from  a  direction  of  100°. 

The  wave  shoaling  derived  from  linear  theory  due  to  both  the 
observed  currents  and  topography  was  calculated  over  the  entire 
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Fig.  tl.  H. s  (m)  from  linear  theory  shoaling  due  to  topography  and  currenis, 
assuming  col  inear  waves  and  currents.  H >  contours  shown  every  0.02  m  from 
0.62  to  0.7  m.  (a)  Incoming  waves  from  CHLV2  at  0. 17  Hz  in  the  wind  direction 
120°.  (b)  0.07  Hz  swell  component  with  direction  100° ,  Hs  contours  shown 
every  0,2  m  from  0.5  to  0.62  m. 


radar  domain  (Fig.  1 1)  for  both  the  wind  wave  and  swell  com¬ 
ponents,  There  was  little  influence  of  topography  for  the  wind 
waves  (0.17  Hz)  because  the  depths  were  >5  m  and  the  waves 
at  this  time  were  quite  small.  There  was.  however,  a  ^25%  in¬ 
crease  of  Hs  due  to  the  surface  current.  The  radar  results  show 
the  same  order  of  effect  localized  to  the  region  of  the  estuarine 
outflow.  The  swell  wave  (0,07  Hz)  was  not  significantly  influ¬ 
enced  by  the  current,  but  was  slightly  enhanced  over  the  shal¬ 
lower  regions  close  to  the  shoreline  and  near  the  shoal  in  the 
center— north  of  the  radar  domain. 

The  winds  during  the  period  from  303.7  to  305,4  were  consis- 
tcntly  from  the  east  or  southeast  but  varied  in  magnitude  from  5 
to  12  m/s.  On  YD  305.33,  the  wind  was  12.3  m/s  in  the  onshore 
direction  at  the  peak  of  the  ebb  tide.  This  caused  the  near-sur- 
face  currents  to  be  aligned  in  the  wind  direction  over  most  of 
the  domain  Fig.  12.  Near  the  mouth,  the  currents  were  much 
weaker  or  in  the  direction  of  the  outgoing  fide.  At  this  time  there 


Fig,  12,  Ha  (meter)  and  current  velocities  at  peak  of  ebb  tide  on  YD  305.33. 
Velocities  determined  by  the  combination  of  rad  i  a  Is  from  each  site.  Velocity 
scale  is  in  lower  right  comer.  Ha  independently  observed  by  HF-S.  Hs  con¬ 
tours  arc  given  every  0.2  m;  only  values  located  <=26  km  from  HF-S  shown. 
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Fig,  13,  £'  component  of  the  ADCP  velocities  during  COPE-3  from  YD  302 
to  YD  307.  (a)  A I .  <b)  A2.  (c)  A3,  Negative  L  denotes  an  onshore  directed  flow. 


was  an  increased  Hs  observed  in  both  radars  near  the  estuary 
mouth  Fig.  1 2.  The  frequency  spectrum  at  CHLV2  revealed  that 
the  wind  waves  were  dominant  at  this  time  [Fig,  1 0(b)]  with  the 
peak  energy  at  0,15  Hz, 

There  was  significant  vertical  current  shear  southeast  of  the 
estuary  mouth  near  the  A  l,  A2,  and  A3  moorings.  The  vertical 
current  profiles  (Fig.  13)  at  these  locations  exhibited  a  surface 
layer  velocity  in  the  wind  direction  and  a  lower  layer  outflow. 
This  was  in  contrast  to  the  typical  outflow  at  the  peak  of  the 
ebb  tide  which  was  normally  concentrated  in  the  upper  5  m  of 
the  water  column.  At  YD  304.7  and  YD  305,3  at  the  peak  of  the 
ebb  the  outflow  occurred  in  the  middle  of  the  water  column.  The 
near-surface  layer  velocities  at  these  times  were  onshore  in  the 
direction  of  the  wind  forcing. 
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IV.  Discussion 

A .  Val Ida  tion  of  HF  Rada  r  Measu  rem  en  is  Wi ih  In  Situ  Senso  rs 

The  current  velocities  at  the  Chesapeake  Bay  mouth  re¬ 
sponded  to  the  interaction  between  the  tidal  currents,  buoyant 
discharge,  wind  forcing,  and  topography.  The  current  response 
to  wind  forcing  was  further  complicated  by  the  disparate  effects 
of  the  wind  over  the  inner-shelf  and  within  Chesapeake  Bay. 
In  particular,  down  welling  favorable  winds  over  the  shelf  may 
produce  a  water  level  setup  within  the  bay  that  drives  a  strong 
discharge  at  the  mouth  [25].  Comparisons  of  HF  radar  obser¬ 
vations  which  were  averaged  over  approximately  l  km2  with 
point  measurements  were  complicated  by  the  high  horizontal 
current  shear.  However,  it  was  necessary  to  make  these  compar¬ 
isons  to  evaluate  the  accuracy  of  the  remote  observations  and  to 
determine  the  spatial  scales  over  which  the  point  measurements 
were  relevant. 

The  agreement  between  the  HF  radar  current  measurements 
and  the  near  surface  bin  of  the  ADCPs  was  excellent  throughout 
the  experiment  at  moorings  A2-A5.  The  HF  radars  and  the 
ADCPs  responded  similarly  to  the  range  of  currents  observed 
(Fig.  4)  with  regression  slopes  close  to  1 . 

The  resolution  of  a  radial  velocity  estimate  was  0.022  m/s 
based  on  the  pulse  repetition  rate  and  the  sampling  duration. 
The  combination  of  the  two  radial  components  to  estimate  the 
vector  velocity  introduced  additional  error  [1 1  ].  This  error  am¬ 
plification  was  less  than  a  factor  of  2  (0.04  m/s),  except  for  the 
North-South  component  at  A4  and  A 5  and  the  Bast- West  com¬ 
ponent  at  Al,  where  it  reached  a  value  of  2.5  (0.05  m/s). 

The  HF  radars  and  ADCPs  were  operating  within  their  ex¬ 
pected  accuracies.  The  rms  differences  for  V  ranged  from  0.068 
to  OT  13  m/s  and,  for  V\  they  ranged  from  0.079  to  0. 101  m/s 
(Table  I).  These  rms  differences  were  of  the  same  order  as  the 
0,078  to  0. 1 2  m/s  observed  by  [26]  and  were  less  than  the  values 
of  0.1 12  to  0.158  m/s  observed  by  [27],  both  of  which  used 
the  same  OSCR  system  as  this  paper.  The  in  situ  measurements 
for  [26]  were  provided  by  moored  instruments  at  depths  of  4  m 
below  the  surface  and,  for  [27],  they  were  at  a  depth  of  9,5  m. 
The  deeper  in  situ  depths  relative  to  the  2-m  ADCP  upper  bin 
used  in  this  paper  should  have  increased  the  contribution  to  the 
differences  by  Stokes  drift  and  wind-induced  shear  [  1 1],  These 
contributions  totaled  0.025  m/s  in  [II]  and  this  provides  an  in¬ 
dication  of  the  strength  of  this  effect  in  the  current  paper. 

While  baroclinicity  can  be  very  important  in  the  region  of  the 
estuarine  outflow,  typical  layer  depths  of  the  buoyant  flow  are 
~5  m  and  its  influence  should  be  much  larger  at  the  inshore 
moorings.  However,  with  the  exception  of  the  0.046-m/s  bias  in 
[/  at  A2  the  statistics  of  the  comparisons  are  similar  at  all  of 
the  ADCPs.  The  biases  at  A3,  A4,  and  A5  were  <0.02  m/s  for 
both  components  and  the  slopes  ranged  from  0.9  to  1 .2  (Table  T). 
These  comparisons  demonstrated  the  accuracy  of  the  HF  radar 
observation  over  much  of  the  study  region. 

Two  multi  frequency  coastal  HF  radars  (MCR)  and  two  Sea- 
Sonde  HF  radars  also  sampled  the  same  general  region  during 
the  COPE-3  experiment.  The  MCR  was  operated  at  frequencies 
of  4.8,  6.8,  13.4,  and  21.8  MHz  [28].  The  SeaSonde  was  oper¬ 
ated  at  a  frequency  of  25  MHz  [29].  The  most  relevant  MCR  fre¬ 
quency  for  comparison  with  the  OSCR  was  21 ,8  MHz.  The  rms 
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differences  between  the  MCR  at  2 1 .8  MHz  and  the  ADCP  moor¬ 
ings  were  similar  to  those  observed  here  (0,053-0.15  m/s).  The 
MCR  at  A2  and  A3  had  larger  rms  differences  than  the  OSCR 
when  compared  with  the  ADCP  while  the  V  component  off¬ 
shore  had  lower  mis  differences.  Because  of  the  relatively  small 
number  of  points  used  in  the  MCR  comparison  (103^103  sam¬ 
ples)  it  is  not  possible  to  draw  any  firm  conclusions  about  the 
generally  small  differences. 

The  rms  differences  between  the  SeaSonde  radial  velocities 
with  antenna  beam  partem  corrections  and  the  ADCP  moorings 
ranged  from  0.105  to  0.133  m/s  [29].  These  values  were  larger 
than  the  vector  component  rms  differences  observed  by  OSCR. 
Only  radial  SeaSonde  comparisons  have  been  published  so  it 
was  not  possible  to  discern  how  well  the  system  as  a  whole  was 
sampling  the  vector  currents.  Clearly,  however,  the  errors  can 
only  increase  when  combining  two  nonorthogonal  radial  mea¬ 
surements.  By  spatially  smoothing  their  results,  the  SeaSonde 
radial  rms  differences  were  reduced  from  0.071  to  0.096  m/s 
which  was  of  the  same  order  as  the  OSCR  vector  component 
differences. 

Both  the  MCR  and  the  SeaSonde  radial  velocities  were  ob¬ 
tained  with  direction-finding  algorithms  corrected  for  measured 
antenna  beam  patterns.  The  OSCR  phased-array  results  were  not 
corrected  for  local  perturbations  to  the  antenna  beam  pattern.  In 
spite  of  this,  there  was  similar  agreement  between  the  three  sys¬ 
tems  and  the  in  situ  currents.  This  demonstrated  that  the  higher 
SNR  obtained  by  the  phased-array  mitigated  the  need  for  an¬ 
tenna  beam  pattern  corrections. 

Inspection  of  the  time  series  at  A2  (Fig.  4)  revealed  that 
the  differences  between  the  HF  radar  and  ADCP  were  not 
uniformly  distributed  relative  to  the  phase  of  the  M2  tide,  The 
largest  differences  consistently  occurred  at  the  peak  positive 
velocity  (maximum  ebb)  and  occasionally  occurred  at  the  peak 
negative  velocity  (maximum  flood).  These  differences  were 
typically  ^0.10  m/s  but  on  several  occasions  were  ^0,20  m/s. 
These  large  differences  at  the  peak  ebb  contributed  significantly 
to  the  0.046-m/s  bias  in  the  U  component  as  there  was  little 
increased  surface  velocity  in  the  V  component  and  only  a 
0.01 8-m/s  bias.  These  differences  were  also  not  observed  at  the 
other  locations. 

Previous  studies  have  demonstrated  that  the  M2  tide  is 
primarily  barotropic  [16].  With  the  exception  of  frictional  influ¬ 
ences  near  the  bottom  it  is,  therefore,  likely  that  local  topography 
may  have  affected  the  comparison  at  this  location.  The  A2 
mooring  was  located  at  the  inshore  edge  of  the  deep  shipping 
channel  leading  into  Chesapeake  Bay  (Fig,  3).  At  this  mooring, 
the  flow  exhibited  significant  spatial  variability  particularly 
during  a  strong  outflow  (Fig,  9).  The  HF  radar  observed  the 
velocity  over  a  footprint  that  was  nominally  1  km  in  horizontal 
extent.  The  cross -she If  scale  of  the  convergence  front  of  an 
estuarine  outflow  is  typically  an  order  of  magnitude  smaller  than 
the  HF  radar  resolution  [30].  Because  of  this  relatively  large 
footprint,  the  HF  radar  observations  may  have  included  current 
velocities  within  the  channel  that  were  not  observed  at  A 2. 

B.  Comparison  of  Wave  Measurements  From  HF-N  and  HF-S 

There  are  expected  differences  in  the  wave  observations  be¬ 
tween  the  point  measurement  from  the  laser  height  gauge  on 
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TABLE  HI 

Mean  and  Variance  of  H*  Observations 


Sensor 

Samples 

Mean  Hs 
[m] 

Var  Hs 
[m!] 

CHLV2 

795 

0,97 

027 

HF-S 

1064 

1.04 

0,32 

HF-N 

898 

1.09 

0.29 

CHLV2  and  the  spatially  averaged  observations  from  the  radars. 
Comparisons  are  further  complicated  by  the  dropouts  at  CHLV2 
at  large  Hs  and  the  spikes  in  the  HF-N  record  at  low  Hs,  The 
challenge  was  to  intercompare  these  different  types  of  observa¬ 
tions  and  assess  whether  observed  differences  in  wave  parame¬ 
ters  were  statistically  significant.  Krogstad  et  al  [3 1  ]  presented 
optimal  regression  methods  to  intercompare  different  types  of 
wave  measurements  that  are  followed  here  for  that  purpose* 

The  most  commonly  applied  statistical  tests  for  intercompar¬ 
ison  of  measurements  are  one-sided  regressions  as  done  for  the 
near-surface  current  measurements.  The  one-sided  regressions 
between  the  radar-derived  estimates  and  CHLV2  revealed  that 
during  the  times  when  there  were  valid  data*  HF-S  had  better 
agreement  than  HF-N  (Table  II).  In  fact*  HF-S  agreed  better 
with  CHLV2  than  with  HF-N. 

If  the  variances  of  the  time  series  are  similar  as  was  the  case 
in  this  paper  (Table  III),  then  it  is  more  appropriate  to  perform 
a  two-sided  regression.  The  two-sided  regressions  through  the 
origin  have  the  same  slope  in  all  three  cases  [Fig.  6(a)-(c)];  how¬ 
ever*  it  was  clear  that  the  comparisons  of  HF-N  with  the  other 
two  measurements  contained  larger  rms  differences  (Table  IV). 
The  statistics  of  the  comparison  between  HF-S  and  CHLV2 
were  within  expected  ranges  based  on  the  variability  of  the  two 
systems*  The  comparison  of  HF-N  to  HF-S  and  CHLV2  re¬ 
vealed  higher  rms  differences  and  a  larger  bias  than  has  been  ob¬ 
served  in  previous  comparisons  between  in  situ  measurements 
[32],  The  biases  from  the  one-sided  regression  involving  HF-N 
were  also  larger  than  those  observed  by  [3 1  ]  and  133]  when  com¬ 
paring  HF  radar  observations  with  in  situ  measurements.  How¬ 
ever,  the  two-sided  regressions  and  ML  estimates  were  similar 
to  those  reported  in  [32]  (Table  IV). 

The  maximum-likelihood  (ML)  approach  as  presented  in  [311 
was  used  to  determine  how  the  sensor  intercomparison  varied 
with  He*  The  percent  relative  error  as  a  function  of  Hs  de¬ 
cayed  for  each  of  the  comparisons  [Fig.  6(d)].  However*  there 
was  much  higher  relative  error  at  low  waveheights  for  the  cases 
involving  HF-N.  This  resulted  from  the  intermittent  large  spikes 
in  the  Hs  as  observed  from  this  station.  For  values  of  Hs  <1 
itl  the  relative  error  of  HF-N  versus  CHLV2  or  HF-S  increased 
rapidly  from  10%  to  >  25%  for  Hs  <0*5  m. 

For  large  waveheights*  the  relative  error  for  the  comparisons 
with  CHLV2  remained  at  >5%  while  the  HF-N  to  HF-S  rel¬ 
ative  error  decreased  rapidly.  This  higher  error  was  likely  the 
result  of  the  dropouts  in  the  laser  measurements  at  larger  wave- 
heights  (Fig*  5).  At  the  peak  of  the  storm  on  YD  293  (Oc¬ 
tober  20*  1997)  the  largest  Hs  recorded  by  CHLV2  was  3*4  m 
while  other  National  Data  Buoy  Center  (NDBC)  wave  measure¬ 
ments  north  (buoy  N44009)  and  south  (DUCN7)  of  the  station 
recorded  values  of  3.9  and  4.0  m,  respectively.  The  NDBC  buoy 


60  ttm  offshore  the  measurement  region  (N440I4)  was  not  oper¬ 
ating  at  the  time  of  the  storm.  Furthermore,  there  was  only  one 
value  above  3.5  m  recorded  at  CHLV2  during  the  entire  year  of 
1997.  These  factors  demonstrate  that  the  CHLV2  laser  height 
sensor  was  not  reliable  for  large  waveheights  (Hs  >3  m). 

The  time  series  (Fig*  5)  of  the  three  colocated  Hs  measure¬ 
ments  shows  that  the  problems  with  HF-N  were  not  systematic. 
Rather*  there  were  frequent  but  intermittent  periods  in  which 
HF-N  overpredicted  Ns *  Inspection  of  the  Doppler  spectra  at 
HF-N  during  the  period  when  there  were  large  differences  be¬ 
tween  HF-N  and  HF-S  and  CHLV2  revealed  large  peaks  at 
frequencies  far  removed  from  the  Bragg  peaks  (Fig.  7)*  These 
spectra  were  not  qualitatively  similar  to  the  normal  spectral  form 
(Fig*  2),  Consequently,  this  and  other  techniques  to  extract  the 
waveheight  from  the  Doppler  spectrum  will  fail.  Because  of  the 
prevalence  of  the  occurrence  of  these  spurious  spectra  in  the 
HF-N  record*  care  must  be  taken  to  eliminate  these  values  when 
interpreting  the  Hs  values  obtained  at  this  site. 

The  noise  in  the  spectra  observed  at  HF-N  was  likely  caused 
by  a  local  or  directional  source  which  was  not  observed  in  HF-S* 
or  from  a  distant  source  that  was  more  strongly  coupled  to  the 
received  signal  at  HF-N  than  HF-S  because  of  the  particular 
antenna  characteristics  at  that  site.  Although  it  was  not  possible 
to  determine  precisely  which  phenomena  was  the  cause  here* 
it  was  likely  that  a  local  source  was  responsible*  HF-N  was 
located  in  close  proximity  (~300  m)  to  a  U.S.  Navy  commu¬ 
nication  facility  that  intermittently  used  HF  frequencies  likely 
leading  to  the  observed  spurious  signals.  For  this  analysis,  we 
will  focus  on  the  more  reliable  HF-S  observations* 

The  comparison  of  the  HF-S-derived  Hs  with  CHLV2  was 
encouraging  for  the  purpose  of  mapping  variability  over  the  out¬ 
flow  plume.  The  26- km  range  from  HF-S  to  CHLV2  was  near 
the  practical  limit  of  wave  observations  from  the  OSCR  system 
[33]*  This  limitation  is  because  the  radar  waveheight  observa¬ 
tions  require  the  extraction  of  the  energy  in  the  second-order 
spectra  (Fig*  2)*  With  increasing  range,  the  second-order  signals 
are  obscured  by  noise  before  the  first-order  peaks  used  to  derive 
the  radial  currents.  The  plume  frontal  location  varied  consider¬ 
ably*  but  often  a  strong  convergence  front  was  present  10-20  km 
offshore  HF-S.  Therefore*  the  measurements  in  this  region  with 
better  SNR  than  the  comparison  point  should  be  quite  reliable* 

This  empirical  approach  was  surprisingly  robust.  The  co¬ 
efficients  required  to  scale  the  linear  relationship  between 
the  second-order  power  and  Hs  were  not  determined  locally* 
Rather,  they  were  derived  from  an  experiment  three  years 
earlier  at  Duck*  NC*  located  ^50  km  south  of  the  present  study 
area*  This  demonstrated  that  the  measurement  was  not  highly 
sensitive  to  local  perturbations  to  the  antenna  beam  pattern* 
This  gives  confidence  that  HF-S-observed  spatial  variability  of 
Hs  was  real. 

C  Wave-Current  Interaction 

In  light  onshore  winds  and  outgoing  tide*  enhanced  wave- 
heights  were  evident  in  the  region  of  the  main  shipping  channel 
(Fig.  9).  When  the  wind  forcing  was  light,  locally  generated 
waves  were  concentrated  at  high  wave  numbers.  The  frequency 
spectra  recorded  at  CHLV2  on  YD  287.78  revealed  a  large  swell 
component  at  a  frequency  0.07  Hz  as  well  as  the  wind  wave  peak 
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TABLE  IV 

Two-Sided  Regression  and  ML  Method  Parameters 
for  Hs  Comparison 


rms  difference 
Hs  [m] 

Two-sided 

regression 

slope 

Bias  of  ML 
method  [m] 

Slope  of  ML 
method 

HF-N  versus 
CHLV2 

0,26 

1.09 

0.08 

1.01 

HF-S  versus 
CHLV2 

0,20 

L09 

0.02 

1.05 

HF-N  versus 
HF-S 

0.27 

1.09 

0,15 

0.96 

at  0.17  Hz  Fig.  10(a).  On  YD  305.33,  the  wind  wave  peak  was 
significantly  larger  than  any  other  component  of  the  spectrum 
[Fig.  10(b)], 

The  wind  waves  were  not  sensitive  to  shoaling  or  refraction 
due  to  topography.  They  did  not  "feel"  the  bottom  in  our 
measurement  domain  because  it  did  not  get  any  shallower  than 
—5  m.  Conversely,  the  influence  of  current  velocity  on  these 
same  waves  was  significant  (Fig.  9),  because  the  fraction  of 
wave  energy  that  was  affected  was  higher  at  low  winds  (Fig.  1 ). 
The  effects  of  nonlocal  swell  could  be  important  in  any  wind 
regime,  but  would  be  more  pronounced  in  light  w'inds.  Linear 
theory  estimates  of  depth  and  current-influenced  wave  shoaling 
reveal  a  peak  associated  with  the  strong  opposing  current  near 
the  estuary  mouth  [Fig.  1 1(a)].  Estimates  of  the  wave  refraction 
(not  shown)  using  ray  theory  showed  no  focusing  of  wave 
energy  within  the  radar  domain.  Had  the  wind  waves  been 
obliquely  incident  to  the  current  shear,  then  more  refractive 
effects  would  have  been  observed. 

The  swell  component  (0.07  Hz)  was  expected  to  be  less  af¬ 
fected  by  the  current  Fig,  I  and  more  affected  by  the  topography 
than  the  wind  wave  component.  Linear  theory  estimates  of  the 
effect  on  the  waveheight  by  both  the  current  and  topography 
predict  a  small  enhancement  of  the  waveheight  over  the  shal¬ 
lower  regions.  However,  because  of  the  reduced  effect  of  current 
on  low-frequency  waves  there  was  no  increased  Hs  associated 
with  the  strong  opposing  current.  There  was  some  indication  in 
the  r:id:ir  nh  tv  at  i  ons  tha  t  th  e  wave  he  i  g  ht  was  i  ncreasi  ng  i  n  t  he 
shallow  region  inshore  of  the  strong  outflow.  This  was  consis¬ 
tent  with  shoaling  of  the  swell  component, 

Hs  observations  from  HF-S  on  YD  305,33  were  higher  off¬ 
shore  the  southeast  portion  of  the  domain,  smaller  in  the  center, 
and  increased  near  the  estuary  mouth.  The  waveheights  offshore 
the  southeast  were  larger  than  elsewhere  in  the  domain  (Fig.  12). 
This  was  in  a  region  where  the  depths  were  greater  than  10  m 
and  the  tidal  currents  were  weak.  It  was  likely  that  the  increased 
waveheights  in  this  region  resulted  from  the  southeasterly  winds 
that  had  been  blowing  earlier  in  the  day,  A  southeast  wind  over 
the  measurement  region  had  unlimited  fetch  and,  consequently, 
swell  could  have  been  propagating  from  this  direction.  Without 
local  directional  spectrum  information  however,  this  could  not 
be  established.  In  the  middle  of  the  domain,  the  waveheights 
were  smallest.  They  then  increased  by  20-25%  in  the  region 
southeast  of  the  estuary  mouth. 

Weak  near-surface  currents  on  YD  305.33  opposed  the  wave 
propagation  in  the  region  of  increased  Hs ,  They  averaged 


20  cm/s  less  than  the  currents  in  the  wind  direction  farther 
offshore.  Since  this  was  the  peak  of  the  ebb  tide,  it  was  likely 
that  the  wind- forced  near-surface  currents  were  interacting  with 
the  normally  barotnopic  [16]  M2  tidal  outflow.  Inspection  of 
the  ADCP  profiles  at  the  A2  and  A3  moorings  revealed  that  the 
vertical  structure  of  the  tidal  flow  changed  during  the  period  of 
strong  onshore  winds  (Fig.  13).  The  offshore  flow  at  the  peak 
of  ebb  tides  both  before  and  after  this  period  was  strongest  in 
the  upper  5  m.  However,  from  YD  304,5  to  305,5,  the  peak 
ebb  tidal  velocities  shifted  deeper  in  the  water  column  and 
there  was  onshore  flow  in  the  near-surface  layer.  Unfortunately, 
salinity  and  temperature  profiles  were  not  available  at  this  time; 
however,  it  is  likely  that  this  was  an  unstable  situation  with 
denser  shelf  water  being  forced  over  the  outgoing  estuarine 
flow.  This  region  then  would  have  been  undergoing  significant 
mixing,  particularly  given  the  relatively  large  waveheights. 

The  waveheights  were  reduced  where  the  surface  layer  was 
moving  in  the  wind  direction  (Fig,  11),  Kirby  and  Chen  [34] 
derived  a  hyperbolieally  decaying  expression  for  the  equiva¬ 
lent  current  velocity  for  a  depth- varying  current  at  a  given  wave 
number.  For  the  wave  numbers  present  at  this  lime  Fig.  10(b), 
their  expression  predicted  little  influence  of  the  observed  cur¬ 
rents  on  the  waves.  Our  observations  were  consistent  with  [34] 
in  that  it  was  only  near  the  estuary  mouth  and  close  to  shore 
where  the  surface  currents  were  opposing  the  wind  direction 
[Figs.  1 1  and  13(c)]  that  enhanced  waveheights  were  observed. 

For  most  users  of  near-real  time  wave  information,  it  is  not 
necessary  to  describe  the  surface  waves  in  great  detail.  Rather, 
significant  waveheights  and  mean  directions  are  most  likely  to 
be  of  primary  interest.  It  is  clear  from  the  results  presented  here 
that  phased-array  radars  can  provide  information  on  the  spatial 
variability  of  This  could  be  of  significant  interest  in  lo¬ 
cations  such  as  the  study  area  where  strong  currents  can  lead 
to  locally  enhanced  Hs .  The  difference  between  an  Hs  of  1,4 
and  2,0  m  over  a  few  kilometers  [Fig.  11(a)]  may  be  signifi¬ 
cant  for  safe  operation  of  small  vessels.  With  the  position  of 
the  outflow  plume  highly  dependent  upon  the  discharge,  tide 
stage,  and  wind  [35],  it  is  not  possible  for  a  few  point  measure¬ 
ments  to  adequately  sample  the  spatial  variability  of  Hs  over 
the  outflow  plume.  Since  the  Hs  estimates  from  the  radars  are 
empirically  based,  it  is  necessaiy,  however,  to  have  at  least  one 
reliable  in  situ  wave  observation  for  calibration  and  verification 
purposes. 

To  establish  a  predictive  capability  for  the  local  wave  con¬ 
ditions,  it  is  necessary  to  implement  linked  shallow- water  cir- 
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culation  and  wave  models.  Such  an  approach  would  be  greatly 
enhanced  by  the  observation  of  the  wave  directional  spectra 
as  in  [33].  Given  the  requirement  for  in  situ  calibration  and 
long-term  observation  over  a  range  of  wind  and  current  condi¬ 
tions  and  the  considerable  expense  associated  with  these  efforts, 
it  is  clear  that  the  necessary  information  to  extract  the  full  direc¬ 
tional  spectrum  from  phased-array  systems  should  be  retained 
whenever  possible. 

V.  Summary 

Surface  currents  and  significant  waveheights  have  been  suc¬ 
cessfully  observed  with  an  HF  radar  system  deployed  near  the 
mouth  of  the  Chesapeake  Bay.  The  remotely  sensed  currents  and 
waves  were  compared  with  in  situ  observations  using  ADCPs 
and  a  laser  height  gauge,  respectively.  The  agreement  between 
the  HF  radar  observed  near-surface  currents  and  the  near-sur¬ 
face  bin  of  the  ADCPs  was  very  good  throughout  the  experi¬ 
ment.  Observed  differences  were  similar  to  or  less  than  previous 
comparisons  between  HF-radar-de rived  velocities  and  in  situ 
point  measurements.  In  particular,  the  vector  component  differ¬ 
ences  were  of  the  same  order  as  beam  pattern  corrected  radial 
differences  from  a  direction-finding  system  deployed  during  the 
same  experiment.  This  demonstrated  that  antenna  beam  pattern 
corrections  were  not  necessary  for  obtaining  reliable  surface 
currents  from  the  phased-array  system. 

Comparison  of  radar  Hs  measurements  with  an  in  situ  laser 
revealed  differences  of  the  same  order  as  in  earlier  studies.  Al¬ 
though  there  were  conditions  under  which  one  of  the  radar  units 
estimated  erroneously  large  Hs,  the  slope  of  the  regression  be¬ 
tween  the  wave  observations  was  only  1 .09.  These  intermittent 
large  values  resulted  from  broad  spurious  peaks  in  the  eeho- 
Doppler  spectrum  which  could  be  removed  a  priori  Given  that 
the  comparisons  were  made  near  the  practical  range  limit  (25 
km)  of  wave  observations  with  the  OSCR,  the  agreement  with 
in  situ  point  observations  provided  confidence  that  the  observed 
spatial  structure  of  Hs  was  valid. 

The  combined  wave  and  current  information  over  this  tidal  ly 
dominated  region  with  complex  topography  revealed  regions 
with  significant  wave-current  interaction.  There  were  several 
occasions  during  the  45  d  of  the  experiment  when  the  incident 
waves  were  propagating  into  an  opposing  current.  The  wave- 
heights  in  low- moderate  winds  increased  by  25%-50%  above 
incident  values,  which  was  comparable  to  linear  theory  esti¬ 
mates,  When  there  was  stronger  wind  forcing,  the  near- surface 
currents  were  directed  w  ith  the  wind,  while  the  tidal  outflow  cur¬ 
rents  below  3  m  were  in  opposition  to  the  wind  direction.  This 
probably  caused  substantial  mixing  between  estuary  and  shelf 
water. 

Real-time  dissemination  of  the  locations  where  waves  are 
enhanced  could  improve  the  safety  of  commercial  and  recre¬ 
ational  maritime  activities.  Longer  term  measurements  and 
additional  comparisons  with  models  for  wave  refraction  and 
diffraction  in  coastal  regions  will  be  necessary  to  quantify  the 
effects  of  opposing  currents  on  the  wave  field  over  a  range  of 
conditions.  Also,  derivation  of  the  full  directional  wave  spec¬ 
trum  from  phased-array  radar  measurements  would  provide 
more  comprehensive  information.  This  could  then  be  incorpo¬ 


rated  into  regional  wave  and  circulation  models  to  develop  a 
predictive  capability. 
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